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Abstract. A novel approach towards the hadronic contributions to the anomalous magnetic moment of the muon is 
presented, namely the Dyson-Schwinger equations of QCD. It has the advantage of being valid for all momentum scales and 
has the potential to address off-shell amplitudes. We present our first results for the pseudoscalar (PS) meson exchange and 
the quark loop contributions. The meson exchange (71°, rj, rj'), aJ; BL ' PS = (84 ± 13) x 10~ n , is commensurate with previous 
calculations, while the quark loop contribution a L B M uarkloo P _ ^Q7±48) x 10 -11 , is strongly enhanced by vertex dressing 
effects in the quark photon vertex. Taken seriously this leads to the estimate of = 116591865.0(96.6) x 10 , giving a 
1.9 <J deviation between theory and experiment. 
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INTRODUCTION 

The anomalous magnetic moment a^ = (g^ — 2)/2 of the 
muon is one of the most precisely determined quantities 
in particle physics. Experimental efforts at Brookhaven 
and theoretical efforts of the past ten years have pinned 
flu down to the 10~ n level, leading to significant devia- 
tions between theory [1] and experiment [2] of ~ 3 a: 

Experiment: 1 16592089.0(63.0) x lO" 11 , (1) 
Theory: 1 16591 790.0(64.6) x 10~ n . (2) 

This discrepancy makes a^ very interesting since it 
might be taken as a hint towards physics beyond the stan- 
dard model (SM). In order to confirm this hypothesis the 
uncertainties of both theory and experiment have to be 
reduced even further. 

The theoretical error is dominated by hadronic contri- 
butions, e.g diagrams that involve QCD beyond pertur- 
bation theory (see [1]). The leading QCD contribution is 
given by the hadronic vacuum polarisation (HVP) inser- 
tion shown in Fig. 1 (left). This diagram also dominates 
the present error of the SM prediction. However, since 
the HVP can be related to experimental data from mea- 
surements of e + e~ — > hadrons, a systematic reduction of 
the uncertainty will come with improved experimental 
input. 

When this is realised, another contribution will dom- 
inate the theoretical error. This is the hadronic light-by- 
light (LBL) scattering shown in Fig. 1 (right). For this 
contribution one must rely entirely on theory since no di- 
rect experimental constraints are available. This contri- 




FIGURE 1. Classes of corrections to the photon-muon ver- 
tex function: (left) Hadronic vacuum polarisation contribution; 
(right) Hadronic light-by-light (LBL) scattering contribution. 



bution involves the hadronic photon four-point function 
which is the central object in the following considera- 
tions. This Green's function has been the subject of inten- 
sive investigation in the past, mainly from the perspec- 
tive of large-A^ e and chiral effective theory. These suggest 
an ordering of diagrams that serve as an approximation 
towards the full function [3]. The diagrams have been 
calculated in the extended Nambu-Jona-Lasinio model 
(ENJL) [4, 5] and the hidden local symmetry model 
(HLS) [6]. Later a refined analysis of the potentially 
leading n° exchange contribution based on ideas of vec- 
tor meson dominance (VMD) has been started [7, 8, 9]. 
Therein, the off-shell behaviour of the % — > yy transition 
form factor was considered. Recently, an analysis of LBL 
within the non local chiral quark model (NL^QM) has 
been presented [10]. 

The hadronic LBL contribution to a^ involves a two- 
fold integration of the four point function (Fig. 1 (right)) 
which makes it a two-scale problem. Thus a separa- 
tion of hard and soft scales is in general not straightfor- 
ward. We circumvent the process of matching an ultra- 



violet (UV) description in terms of perturbative quarks 
and gluons and effective mesonic degrees of freedom in 
the infrared (IR) by choosing a description that is en- 
tirely based on the fundamental degrees of freedom of 
QCD, quarks and gluons. In this non-perturbative ap- 
proach, bound-states arise dynamically. We rely on the 
Dyson-Schwinger equations (DSE) which provide a non- 
perturbative means to obtain full one particle-irreducible 
(1PI) Green's functions [12]. In particular we work in a 
truncation scheme defined in Ref. [13] (see the next sec- 
tion). Further details can be found in [14, 15] and exten- 
sions beyond this truncation can be found in [16]. 

In the present model the photon four point function 
may be defined as a resummation of an infinite subset 
of the leading order l/N c diagrams. In particular only 
planar diagrams without internal quark lines are con- 
sidered. Additionally, the Yang-Mills sector of QCD is 
truncated to two-point functions so that only 'rainbow- 
ladder' gluon insertions are taken into account. This 
leaves us with the expansion shown in Eq. (3) 




where all quarks and quark-photon vertices are fully 
dressed. It is well known that the second diagram of this 
expansion includes various resonances and in particular 
pseudo-scalar mesons. This gives rise to the pole approx- 
imation shown in Eq. (4) 




which becomes exact on the respective meson mass shell. 
This leads us to a picture that includes meson exchange 
terms in analogy to the 1 /N c picture. While an evaluation 
of Eq. (3) is currently underway we present here results 
for Eq. (4) as a starting point of our investigations. 



PROPAGATOR AND VERTICES WITHIN 
THE DSE/BSE-APPROACH 

In the following we summarize our calculational scheme 
for LBL; details will be given elsewhere [14, 15]. The 
three objects that we must calculate are the fully dressed 
quark propagator (Fig. 2), the Bethe-Salpeter-amplitude 
(BSA, Fig. 3) for pseudo-scalar (PS) mesons and the 
quark-photon-vertex (Fig. 4). These are defined via self- 
consistent integral equations that require numerical solu- 
tion. 

The gluon (D liV {k 2 )) and quark-gluon vertex are 
needed to solve the quark DSE of Fig. 2. In our trun- 



cation the vertex is considered to be projected onto one 
component r v (p,q) := T YM (k 2 )y v giving the interaction 

•jwssosbbbO = Y^D^ v (k 2 )T YM (k 2 )y v . (5) 

The dressing Z(k 2 ) of the landau gauge gluon D^ v (k 2 ) = 
(S^v — kpkv /k 2 )Z(k 2 )/k 2 is combined with the vertex 
dressing F YM (k 2 ) and modeled by a single dressing func- 
tion that approaches one loop perturbative QCD (pQCD) 
in the UV and contains a strongly enhanced IR part that 
induces dynamical chiral symmetry breaking (D^SB) 
[13]. The resultant quark contains pQCD above a few 
GeV and develops an enhanced mass function of typi- 
cally a few hundred MeV in the IR. In this way perturba- 
tion theory is unified with the constituent quark picture 
with the quark defined at all momentum scales. 

Once the quark-DSE is solved, the quark, together 
with (5), forms an integral part of the remaining two 
equations for the three-point functions (Figs. 3 and 4). 
This truncation respects the chiral properties of QCD 
such that low energy theorems like the Gell-Mann- 
Oakes-Renner relation [17] are fulfilled. In this way ob- 
servables such as bound-state masses, leptonic decay 
constant and form factors in the pseudo-Goldstone octet 
are reproduced on the percent level while the accuracy is 
five to ten percent for vector mesons [13]. These channels 
are the most important ones in the present calculation. 

These building blocks now allow for a description 
of QCD bound-states as solutions of the BSE (Fig. 3), 
together with electromagnetic properties which are de- 
pendent on the quark-photon- vertex (Fig. 4) that couples 
QCD to QED. We emphasize that the self-consistent 
quark-photon-vertex contains a dynamically generated 




FIGURE 2. Dyson-Schwinger equation for the quark prop- 
agator. 

FIGURE 3. Bethe-Salpeter equation for quark-antiquark 
bound states. 




FIGURE 4. Inhomogeneous BS equation for the quark- 
photon vertex. 



vector-meson bound state pole in the time-like region 
thus encapsulating the ideas of vector meson dominance 
(VMD). More details can be found in [18]. Consider as 
an example the PS — > yy transition form factor, shown 
in impulse approximation in Fig. 5. This description ful- 




FIGURE 5. The PS -> yy transition form factor. 

fills the most important constraints for the case of the 
pion, such the charge conservation in the correct asymp- 
totic limits for high photon virtualities. More information 
about the properties of this form factor can be found in 
[19]. In the present case, however, the pion is far away 
from its mass-shell and thus an off-shell prescription is 
necessary (see [15] for details). 

With the quark, quark-photon vertex and PS yy 
specified we now have all of the ingredients necessary 
to calculate the photon four-point function and its expan- 
sion depicted in Eq. (4). 



RESULTS AND DISCUSSION 

The n, 77, 77' exchange contributions amount 



to a 



LBL;tt 1 







(57.5 ± 6.9) x lO" 11 , c% 



LBL;r) 



(15.8±3.5) x 10" 11 and Off Uri ' = (11.0±2.4) x 10" 11 
leading to 



3 LBL;PS = 3 ± 12 x 10 -11 



(6) 



The errors include numerical uncertainties as well as esti- 
mates of the model dependence. The latter is determined 
by consideration of meson observables in the respective 
channels. The result given in Eq. (6) compares well to 
previously obtained values [6, 4, 7]. 

Due to numerical complexity, our results for the quark 
loop contribution do not yet include the full quark- 
photon vertex (Fig. 3). Instead we quote here results 
which include (a) bare vertices and (b) the leading part 
of the Ball-Chiu construction (1BC) [21]: 



LBL;quarkloop (bare vertex) 



t 



BL;quarkloop(lBC) 



( 61 ±2) X 10" 11 

(107±2) x 10" 11 



(7) 



The result for bare vertices in Eq. (7) is comparable 
to that found in constituent quark models. The 1BC 
vertex strongly enhances the contribution to the muon 
anomaly. This is in stark contrast to what has been found 
in other models [4, 6]. The difference of the two results 
in Eq. (7) determines the leading error quoted above. 
The inclusion of all non-perturbative covariant of the 



quark-photon vertex is desirable, but is an extremely 
elaborate calculation and deferred to future work. At 
face value, this then leads to a revised estimate of the 
total a M = 116591865.0(96.6) x 10" 11 , which reduces 
the difference between theory and experiment to ~ 1 .9 a. 
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